Switching current noise and relaxation of ferroelectric domains 
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We simulate field-induced nucleation and switching of domains in a three-dimensional model of fer- 
roelectrics with quenched disorder and varying domain sizes. We study (1) bursts of the switching 
current at slow driving along the hysteresis loop (electrical Barkhausen noise) and (2) the polariza- 
tion reversal when a strong electric field was applied and back-switching after the field was removed. 
We show how these processes are related to the underlying structure of domain walls, which in turn 
is controlled by the pinning at quenched local electric fields. When the depolarization fields of bound 
charges are properly screened we find that the fractal switching current noise may appear with two 
distinct universal behaviors. The critical depinning of plane domain walls determines the universal- 
ity class in the case of weak random fields, whereas for large randomness the massive nucleation of 
domains in the bulk leads to different scaling properties. In both cases the scaling exponents decay 
logarithmically when the driving frequency is increased. The polarization reverses in the applied 
field as a power-law, while its relaxation in zero field is a stretch exponential function of time. The 
stretching exponent depends on the strength of pinning. The results may be applicable for uniaxial 
relaxor ferroelectrics, such as doped SBN:Ce. 
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I. INTRODUCTION 

Polarization reversal and properties of the accompany- 
ing switching current are relevant for various applications 
of ferroelectric materials. In particular, it is a key pro- 
cess that takes place in electronic devices based on ferro- 
electrics thin films, where a reproducible switching pro- 
cess is required. Apart from the extensive experimental 
and theoretical work (TJ-j^], a complete understanding of 
the domain-wall motion and predictability of the switch- 
ing processes remains elusive in real ferroelectrics, espe- 
cially when several physical parameters, determining ei- 
ther structure or dynamics, are varied. Domain switching 
and relaxation (back-switching) in ferroelectrics can be 
studied by hysteresis loop, switching current, and switch- 
ing current bursts (electrical Barkhausen noise). 

The development of domain structure and motion of 
domain walls during switching can be regarded in the 
scope of nonlinear dynamic processes occurring in driven 
disordered systems, which are studied extensively in the 
past decade. Particularly, the measurements on fer- 
romagnets (a short review is given in W) show that 
Barkhausen noise exhibits interesting fractal properties, 
that can be characterized by a full set of scaling expo- 
nents of the sizes, durations and energies of the bursts 
and of the power spectrum. An example with a com- 
plete set of measured scaling exponents and correspond- 
ing fractal dimensions with derivation of the associated 
scaling relations can be found in Ref. H. A possible 
grouping of the scaling exponents into different univer- 
sality classes of scaling behavior offers a new and simple 
link between the measured Barkhausen noise and the do- 
main structure. For example, in a recent study |6| such a 
relation was established between the universal scaling ex- 



ponents of acoustic emission noise and the packing prop- 
erties near the first-order structural transformation. 

The electrical Barkhausen noise has been 
served experimentally in BaTiC>3 flJI]^ 
Gd 2 (Mo0 4 ) 3 and LiTa0 3 g§ single crystals 
in 9.4/65/35(Pb,La)(Zr,Ti)0 3 (PLZT) ceramics 
and PbMg 1/3 Nb2/3C>3 (PMN) relaxor ferroelectric 
Physically, the polarization reversal and domain wall 
motions in ferroelectrics is accompanied with the local 
change of sign of bound charges, which influences values 
of the electric field in the sample. The depolarization 
field of bound charges can be partially screened by an 
adequate addition of charges on the surface and by other 
methods (see [[| and references therein). Hence at a 
given point x in the sample at time t un-compensated 
depolarization field is given by H 



5E dep (x,t) 



E dep (x,t) -^E s 



(1) 



where E sreen (x, t) is a total field due to screening 
(surface and bulk) charges at that site. Incomplete 
screening of the depolarization field tends to change 
the electric field inside the sample, thus influencing the 
course of switching process itself. For instance, there 
is a finite probability that inverse jumps of polariza- 
tion (negative Barkhausen pulses) occur, that slow 
down the switching. Recently a series of sophisticated 
measurements were done [p|JT^| that reveal the effects 
of the depolarization field on a created plane domain 
wall in Gd2(MoC>4)3 and on a multi-domain structure 
in LiTaOs. The obtained distributions of lengths, sizes 
and rest time of Barkhausen jumps obey a power-law be- 
havior with an exponential cut off, both for switching 
and back-switching processes. Using a suitable experi- 
mental set-up the power-spectrum in PMN relaxors was 
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measured, which was found to decay with frequency lu as 

Recently studied uniaxial relaxor ferroelectrics 
Sr .6iBa .39Nb 2 O 6 (SBN61) and Pbo.eiBao.sgNbzOg are 
regarded as systems with the Ising symmetry of the order 
parameter, having the tetragonal structure in the low- 
temperature phase with the spontaneous polarization 
along c-axis HQ. In solid solutions SBN61 : X dop- 
ing with X = Cer+ , Cr 3+ , Co 2+ reduces the transition 
temperature and induces local electric random fields, 
which give rise to local polar clusters Recently a 

comprehensive series of measurements in SBN:Ce and 
similar systems [|l3|-|l5| gave a strong support to the 
presence of random electric fields in this class of relaxor 
ferroelectrics. In general, more common sources of uncor- 
rected random fields in ferroelectrics are neutral defects 
that can locally break the polar symmetry, as for instance 
in TGS (see detailed discussion in Ref. |fT6fl , page 149). 
[Note that, in contrast to ferroelectrics, in ferromagnets 
local random fields are not allowed by the symmetry, but 
are theoretically constructed as being induced by coarse 
graining due to other types of disorder (see also [0)-] 

In this work we simulate temporal behavior of polar- 
ization driven by the electric field in three-dimensional 
systems of Ising spins subject to quenched local random 
fields and positive nearst-neighbor interaction, which is 
compatible with the occurrence of a ferroelectric long- 
range order. We consider two types of driving. First, by 
stepwise increase (decrease) of the external field along 
a hysteresis loop (see below) we study bursts in the 
switching current j(t) = dP(t)/dt (electrical Barkhausen 
noise). We discuss effects of incomplete screening of de- 
polarization fields on the coarse of switching and hystere- 
sis loop properties. When the screening is complete, we 
study how the statistical properties of the noise depend 
on the strength of pinning and on the frequency with 
which the driving field is changed. Second, by imposing 
a high opposite value of the external field in the homoge- 
neous state we let the polarization to reach its saturation 
at that field, and then remove the field. We determine 
reversal and spontaneous back-switching of the polariza- 
tion with time at different strength of disorder. This 
type of switching is more familiar in the experiments on 
ferroelectrics [ |15[ . 

Apart from allowing negative Barkhausen pulses as de- 
scribed below, the present study is complementary com- 
pared to earlier simulations of Barkhausen pulses in 3- 
dimensional Ising systems in the presence of disorder 
JT8|-pH in the following sense: In the earlier work |l8| we 
considered only strong disorder that is compatible with 
nucleation of small domains under driving field, as for in- 
stance in relaxor ceramics, and analyzed the signal from 
only central part of the hysteresis loop, as it is usually 
done on experiments. In the present work the study is 
extended to entire range of disorder including the weak 
pinning region where extended domain walls may occur 



due to lack of pinning centers p2|. In this case a po- 
tential depinning of domain walls at a critical driving 
field contributes significantly to the scaling behavior of 
Barkhausen pulses, as recent theoretical results on Bethe 
lattice suggest Q. Therefore, we include integration 
over entire hysteresis loop in order to incorporate the 
contribution of critical depinning to the universal scal- 
ing properties of the induced electrical noise in the case 
of weak pinning |24| . Furthermore, here we employ the 
simulation algorithm at finite driving rates, that mimics 
closely the experimental situation with stepwise increase 
(decrease) of the electric field. We examine in detail the 
effects of driving rates in both strong and weak pinning 
case by varying height of the steps in the external field. 
The earlier simulations in Refs. 19 2(| were done in the 
limit of zero driving rate, which appears physically in- 
accessible to real experiments. Whereas an entirely dif- 
ferent aspect concentrating on the dynamic scaling dur- 
ing the growth of individual avalanches was considered 
in It should be stressed that in this paper we im- 

plement strictly quenched disorder, as described in Sec- 
tion II. Note that some speed algorithms, as for instance 
used in pO| , pl| , allow certain dynamic variations of pin- 
ning while the avalanche grows, which can be relevant 
to systems with annealed rather than quenched disorder. 
Theoretically, having unfixed pinning centers during the 
avalanche growth may alter the long-range correlations 
in the avalanche-triggering fields and in distances be- 
tween initial points of avalanches, which were shown to 
be closely associated with the appearance of the fractal 
noise in two-dimensional systems p5[ . 

The paper is organized as follows. In Sec. II we intro- 
duce the model and discuss all relevant parameters that 
influence the hysteresis loop properties. In Sec. Ill we 
study in detail scaling behavior of the switching current 
noise under varying disorder and driving conditions. Sec. 
IV contains the results of switching in a high field and 
spontaneous back switching of polarization when the pin- 
ning is varied. The paper contains a short summary of 
the results and discussion in Sec. V. 



II. MODEL AND RELEVANT PARAMETERS 

We consider a model with Ising-type spins situated on 
three-dimensional simple cubic lattice with positive in- 
teractions between neighbor spins and local random field 
at each lattice site p7[ : 

U=~Y J J xy S(x)S(y) - KS(x) S(x) . (2) 

x,y x x 

We assume that the random exchange interactions J x ,y 
are distributed around a positive value J with a narrow 
distribution. The random fields h x are given by Gaussian 
distribution with zero mean and the width / (measured 



2 



in units of J 2 ). As mentioned above, the source of lo- 
cal random fields can be neutral defects that break po- 
lar symmetry [|l6| or charge disorder due to doping H, 
which cause occurrence of polar regions in relaxor ferro- 
electrics. In addition, a random part of the interaction 
between these polar regions contributes to the observed 
glassy behavior of the dielectric response [^6| . Since the 
random fields locally break the symmetry of the order 
parameter, we believe that the dominant disorder effects 
on the dynamics that we consider are due to random 
fields. Therefore, here we neglect the randomness in the 
spin-spin interaction. 

The spin dynamics that we are considering is field- 
assisted and takes part deeply in the ordered phase. 
Hence, we neglect possible temperature effects, both for 
the reasons of keeping the number of parameters finite 
and assuming that they play a secondary role in these 
processes. The dynamics consists of spin alignment along 
a locally dominant field. Thus, when the local field 
at site x and time t, hi oc (x,t) exceeds zero, the spin 
S(x, t) at that site flips and thus aligns along the local 
field. The local field hi oc (x,t) consists of the interac- 
tion and pinning part hi p (x,t) = J2 y JxyS(y,t) + h(x) 
and the driving field E(x,t). The driving field itself is 
given by E(x,t) = E ext {t) - 8E dep (x,t), where E ext (t) 
is the external field and SEd e p(x,t) a non-compensated 
part of the depolarization field at that site [[|. We 
employ a slow stepwise increase of the external field 
E ex t{ti) = —E sat + tiAE at discrete time intervals {ti}, 
while assuming that a non-compensated part of the depo- 
larization field at site x may assist a probabilistic (prob- 
ability parameter b) back-flip of an already aligned spin 
at that site. It should be stressed that the parameter b 
represents the probability that at a site x in time t the 
depolarization field opposes the effect of external filed. 
(The situation where the depolarization field is parallel 
to the external field is less interesting within the present 
numeric implementation of the dynamics where driving 
rate is not fixed, as described below.) Having in mind 
that the depolarization field at a given site depends on 
the actual domain structure and its kinetics, the uniform 
distribution of the stochastic variable b can be regarded 
only as a first step towards modeling of this complex phe- 
nomena. 

Implementation of the dynamics is as follows: 
Quenched random fields of Gaussian distribution (in dou- 
ble precision) are generated and stored. They are kept 
constant until entire hysteresis loop is completed by the 
driving field. The external field is increased by a small 
amount AE and the set of local fields {h loc (x, t)} com- 
puted at each lattice site and stored. The system is up- 
dated in parallel, i.e., with respect to stored set of lo- 
cal fields, and the spin flipped when h w (x, t) + E ext (x, t) 
exceeds zero by a small amount I0~ 10 . When it is as- 
sumed that the depolarization field is not entirely com- 
pensated, the spin flip is executed with a reduced prob- 



ability I — b. This completes one time step of avalanche 
evolution. Then a new set of local fields at affected sites 
(neighbors of just flipped spins) is computed and stored. 
The conditions for spin flips at those sites examined and 
spins flipped according to the above rule. The process 
continues until no more spins are found which satisfy 
the conditions for flip. This completes one Barkhausen 
avalanche. Then the external field is increased again by 
the amount AE. When the whole loop is completed, we 
set another distribution of quenched random fields and 
repeat the driving. The results are averaged over the 
entire ensemble (up to 200) of different random-field re- 
alizations. 

In this model the following relevant parameters ap- 
pear: Width of the random-field distribution /, that de- 
termines strength of pinning and/or domain size; jump 
AE of the external field, which defines the driving rate; 
probability of back-flips b, which is related to the un- 
screened part of the depolarization field. In practice, 
usually an additional parameter is due to the section of 
the hysteresis loop at which the signal is monitored, since 
the properties of the signal vary along the loop. As men- 
tioned above, here we integrate the statistics over the 
entire hysteresis loop in order to capture the contribu- 
tion of the depinning of large domain walls to the scaling 
behavior of noise by weak pinning. 

The hysteresis loop dependence on pinning by random 
fields and on driving rates is shown in Fig. I. By varying 
the strength of pinning / the form of the hysteresis loop 
changes (Fig. la) from the rectangular loop at low dis- 
ordered 3-d sample to the slim loop, when the pinning is 
strong (large /). In the case of low disorder, i.e., at weak 
pinning of domain walls, large domains may form due to 
lack of pinning centers. In real ferroelectrics occurrence 
of large domains are usually accompanied by extended 
plane domain walls [fjHJ^]. On the other side, when the 
pinning is strong, such as in the case of pinning at grain 
boundaries of small grain ceramics |f(it| , many small do- 
mains can be nucleated at spatially distributed centers 
with weak pinning, and their walls can propagate along 
few active sites with weakest pinning, as the driving field 
increases. This makes effectively low-dimensional walls 
moving in a three-dimensional sample. In addition, a 
massive nucleation of domains in the central part of hys- 
teresis imposes a spatial restriction on the growth of later 
domains. Hence we expect a different scaling behavior 
of sizes and durations of induced switched domains in 
the case of weak and strong disorder. Our simulations 
confirm this phenomenological picture. In ferroelectrics 
a continuous variations of the size of domains can be 
achieved by increasing temperature or by applying stress, 
as for instance described in 8/65/35 PLZT ceramics [ p7| . 
It was suggested in Ref. |ll| that the appropriate illu- 
mination of samples can notably reduce the strength of 
pinning in doped relaxor SBN61 : Ce 3+ ferroelectrics. 

We would like to stress that in the algorithm keeping 
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the external field constant during the avalanche evolution 
and increasing it only after the avalanche has stopped 
does not allow to study rest periods between avalanches. 
However, this type of field update leads to a better reso- 
lution of the noise, while it changes the time dependence 
of the driving field (cf. Fig. fc). The driving rate can be 
better defined as the average derivative r =< dE(t)/dt > 
along the loop. In the simulations we find that when back 
flips occur with a finite probability b > they effectively 
reduce increase of the driving field and prolong the du- 
ration of a jump, as well as inducing inverse jumps, as 
shown in Fig. 1 c,d. 

Assuming that the back flips occur very rarely (b — ► 0), 
we obtain the switching current noise which is shown in 
Fig. 2. The train of bursts of the switching current ex- 
hibits a long-range correlations — Fourier spectrum j(uj) 
of the signal is shown in Fig. 2 (top). It is related to 
the power-spectrum via S(uj) = \j(ui)\ 2 , that can be 
directly measured (see fill). Hence, the exponent of 
the power spectrum in the present model is given by 
2 x <fi = 1.48 ± 0.02 for the case of moderate pinning 
/ = 3.6. The power-law behavior of the power spectrum 
suggests that the switching current noise exhibits fractal 
character when certain conditions are met, in particular, 
when the depolarization fields are successfully screened. 



III. FRACTAL PROPERTIES OF THE 
SWITCHING-CURRENT NOISE 

The fractal noise of the switching current have a range 
of universal scaling properties, that can be determined 
more clearly if the field inside the sample is fully con- 
trolled. Therefore, in this section we put b = in order 
to study in detail these scaling properties and to deter- 
mine how they depend on two relevant parameters — the 
domain size and the driving rate. 



A. Varying strength of pinning 

First we fix the driving rate to a small value AE/J = 
1 x 10 -3 and consider how the domain size or strength of 
pinning influence the properties of the switching current. 
In Fig. 3 we show the cumulative distributions of size 
(area) s of switched domains and duration t of switch- 
ing following a single field update, for several values of 
pinning strength /. As the Fig. 3 shows, both switched 
areas and durations of their switching show a wide range 
of values that are described by power-law distributions 
with finite cut-offs (X stands for X = s,t): 



P(X, L) = X 



-(rx-l) 



V(XL- 



(3) 



The scaling exponents r s — 1 and r t — I determine 
the respective slopes of the cumulative distributions of 



switched area and duration of switching (see Fig. 3). The 
fractal dimension D s of size of switched domains and the 
dynamic exponent D t = z in Eq. (|^) determine how the 
corresponding cut-offs scale with the characteristic length 
(here the system size L), Xo ~ L Dx . 

Both the slopes Tx and the cut-offs, related to Dx, 
of the distributions depend on the strength of pinning. 
First, for weak pinning (small values of /, bottom curve 
in Fig. 3) we notice that extended domains, bounded by 
presumably plane domain walls, may occur. By increas- 
ing the field these domain walls move through a weakly 
random medium and eventually depin at a critical field 
Hd(f). Depinning of large domain walls leads to a peak 
in the histogram at large values of s, or a plateau in the 
cumulative distribution as shown in Fig. 3. Before the de- 
pinning is reached a series of small avalanches occurs that 
are power-law distributed. The plateau disappears at a 
critical value of pinning strength f c (r,L) w 2.4, where 
a pure power-law distribution with virtually infinite cut- 
off (controlled only by the system size) occurs. Here the 
measured slopes give t s = 1.74 and r t = 2.14 for the fixed 
driving rate AE/J = 1 x 10~ 3 (see Fig. 3). 

At strong pinning (large /) many small domains are 
being nucleated at isolated points in the bulk and grow 
slowly along a few active sites per time step. The 
maximum sizes of domains are finite and decrease with 
stronger pinning (cf. top two curves in Fig. 3). Here 
the measured slopes give t s — 1.56 and r t = 1.94 for 
size and duration, respectively, suggesting another uni- 
versality class of the scaling behavior, compared to the 
case of weak pinning. The transition to the strong pin- 
ning regime for / > f c is marked by disappearance of the 
critical depinning of large domain walls at all fields and 
domination of the nucleation processes. Numerical data 
in Fig. 3 suggest that the critical value f c is close to 2.4 
for the used system size and driving rate. 

Qualitatively similar behavior with two universality 
classes but with different exponents, was found in the 
case of two-dimensional system where an extended one- 
dimensional domain wall was initially prepared |p8| . For 
comparison, in Table 1 we summarize exponents both for 
3-dimensional and for 2-dimensional case. (Notice that 
the listed exponents for three-dimensional system are ob- 
tained for a small finite driving rate AE/J = 4x 10~ 4 , 
whereas the simulations at zero driving rate were em- 
ployed in two dimensions.) The results for 2-dimensional 
case are expected to be relevant for domains on the sur- 
faces of the bulk samples (see Ref. |14 ) and for thin films 
in which the width of the domain wall exceed the film 
thickness. An interesting physical realization of varying 
strength of pinning was suggested recently pjJ in ferro- 
magnetic C o /C o film by varying temperature through 
the transition temperature at which the substrate or- 
ders an£i-ferromagnetically. As already mentioned, in 
ferroelectrics the local random fields are physical [|l3] . 
Therefore varying the strength of pinning due to random 
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fields in ferroelectric thin films and bulk materials can 
be achieved in a more direct way, offering diverse pos- 
sibilities to detect the change in the scaling behavior of 
the switching current noise and critical properties on the 
hysteresis. 

B. Varying driving rate 

In the following we examine how the scaling properties 
of the switching current noise change when the driving 
rate is increased. Very low frequencies of the driving field 
are usually not accessible on experiments, both for phys- 
ical and technical reasons. A finite driving rate, which is 
defined as r =< dE(t)/dt >, can be approximated with 
AE/J in the case of stepwise increase of the external field 
in the absence of back flips 6 = 0. For finite rate r > 
the switching may start at different parts of the sample 
and the measured signal is a superposition of these sep- 
arate events. In contrast to the mathematical limit of 
infinitely slow driving, by finite field jump AE the do- 
main wall may get enough energy to overcome several 
pinning centers at the same time. Thus either different 
small domains are initiated simultaneously, or different 
parts of a large domain wall advance in parallel. In both 
cases the statistical properties of the signal are changed. 
Merging of many individual jumps makes the occurrence 
of large events more probable. 

For a range of values of driving rates we still find a 
power-law behavior with reduced values of the scaling ex- 
ponents and enlarged cut-offs. In Fig. 4 we show how the 
distributions for sizes and durations of pulses vary when 
the driving rate AE/J is increased, while the strength of 
pinning was kept constant in the strong pinning regime. 
In Fig. 5 the same type of analysis was done for the case 
of weak pinning. In the strong pinning regime a dramatic 
increase of the cut-offs in the distribution of switched area 
occurs due to merging of many small signals. Whereas, in 
the weak pinning the relative enlargement of the cut-off 
is smaller, that can be understood as an effect of parallel 
advancement of different sections of the same extended 
wall, rather than merging of different domains. Both in 
strong and weak pinning case the increased incidence of 
large events leads to decrease of the scaling exponents. 
We find logarithmic dependence of the exponents for the 
size and duration of events on AE/J both in weak and 
strong pinning regimes. The results are shown in Fig. 6 
and the fit lines are given by the general form 

t x = A x - B x \n{AE/J) . (4) 

For the strong pinning both area and duration exponents 
decrease equally, i.e., B s w B t ~ 0.08, with the con- 
stants A s « 1.0 and A t ~ 1.4. In the weak pinning, 
however, the duration exponent decreases with nearly 
double rate compared to the size exponent. We find 



B t « 2B S w 0.11, (see Fig. 6) and A t w 1.37, A s w 1.30 
within numerical error bars. 

It is important to note that the measuring conditions 
depend on the experimental set up and they can vary 
from one experiment to the other. The effective driv- 
ing rate depends on the steps in the external field and 
the degree of compensation of depolarization fields in- 
side the sample. In real experiments the driving rate 
is always finite, resulting in the scaling exponents that 
are always smaller (and the fractal dimensions larger) 
compared to those computed in the theoretical limit of 
zero driving rate r — > 0. Therefore, in order to compare 
measured and predicted scaling exponents it is advisable 
to make several measurements at different driving rates 
and find out if the results, for instance for the area expo- 
nent r s , belong to a line in the range 1.66 — 1.40, which 
is compatible with occurrence of small domains, or in 
the range 1.78 — 1.60, suggesting that the sample has ex- 
tended plane domain walls. When the driving rate is too 
large (cf. top line on both panels in Fig. 5), the scal- 
ing behavior of the distributions is lost. The cumulative 
probability distribution shows curvature both due to the 
lack of small events and the excess of large events. 



IV. POLARIZATION REVERSAL AND 
BACK-SWITCHING 

In the previous section we have shown how different 
kinds of domain walls in the sample determine the scal- 
ing properties of switching current noise, reflecting the 
nature of domain-wall motion when the electric field is 
slowly ramped. Here we consider another type of driv- 
ing. First, in the initially polarized state (all spins down) 
we impose a large opposite field E s and keep it until po- 
larization reverses and reaches saturation in that field. 
Then we switch off the field to zero and monitor how 
the polarization decays. Both processes depend on the 
type of domains in the sample, and thus, on the strength 
of pinning. The results for polarization as a function of 
time are shown in Fig. 7 for various strengths of pinning 
/• 

For weak pinning spontaneous back switching of the 
polarization is exponentially fast, reaching the value P r 
of remanent polarization. Whereas, the relaxation be- 
comes increasingly slow as the strength of pinning / in- 
creases. The limiting values of P r are lower, which cor- 
respond to increasingly slimmer loops for strong pinning 
(cf. form of the hysteresis loop in Fig. la). In general, 
polarization decay in Fig. 7 can be approximated with a 
stretch exponential law P{t) = ^4 + £> exp[— [x/t)"], with 
a = a(f ) decreasing with the pinning strength. We have 
a wl, 0.92 and 0.80 for / = 2.4, 3.4, and 4.4, respec- 
tively. Qualitatively similar back-switching was recently 
observed in doped SBN : Ce 3+ relaxor ferroelectric with 
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two different stretching exponents corresponding respec- 
tively to illuminated and non- illuminated sample |l5| . 

The polarization reversal in the strong constant field 
shows slow increase in the presence of small domains, 
but with another type of time dependence (see Fig. 7). 
The positive part of the polarization curve is fitted by the 
power-law increase P(t) = 1 — C(t — a) - ", where C « 0.8, 
a is the respective shift (see Fig. 7) and the exponent v k 
2, 1.75, and 1.68, for / = 2.4, 3.4, and 4.4, respectively. 
Again, the characteristic exponent v decreases when the 
pinning becomes stronger, thus contributing to the slow 
dynamics but within a power-law, rather than a stretch- 
exponential dependence. 



V. CONCLUSIONS AND DISCUSSION 

We have shown that, in addition to the measurements 
of the hysteresis loop properties and smooth component 
of the switching current, the domain structure in a ferro- 
electric sample can be indirectly determined by the non- 
invasive dynamic measurements of the switching current 
noise and back-switching processes. Fractal character of 
the dynamic response indicates that a nontrivial domain 
structure can occur under slow held driving. The mea- 
sured properties in these processes can be expressed in 
terms of universal exponents, which turn to be related in 
a unique way to the underlying structure and motion of 
domain walls. We have shown how the universal features 
of the switching current noise are quantified through two 
distinct sets of scaling exponents. The universal scaling 
behavior is associated with a dominant scenario of the 
domain wall dynamics: either the motion and depinning 
of large domain walls, in the case of weak pinning, or 
bulk nucleation and growth of many small 3-dimensional 
domains, in the case of strong pinning. The pinning in 
our model is provided by the local electric random fields. 

Similarly, the polarization reversal in a strong applied 
field, and back-switching of the domains when the field is 
set to zero, are quick when large domains occurs, whereas 
they are slow when many small domains are nucleated. 
In the latter case the polarization relaxes according to a 
stretch-exponential law (the stretching exponent a < 1) 
whereas it reverses with a fractal power of time (the 
characteristic exponent v < 2) when a large opposite 
field is switched on. In the case of weak pinning, on 
the contrary, the relaxation is exponential (cr « 1) and 
the reversal algebraic in time with v = 2. Dependence 
of the stretching exponent cr(f) on the pinning strength 
is in qualitative agreement with measurements of back- 
switching in doped SBN : Ce 3+ relaxor ferroelectrics [[l5| , 
again suggesting that the local electric fields contribute 
significantly to slow the dynamics in this class of relaxors 
p5[ . (Note that the characteristic time appears unreal- 
istically small in our lattice model.) We are not aware 



of measurements of temporal behavior of the polariza- 
tion reversal and of the switching current noise in these 
systems. 

The observed scaling behavior of the switching current 
depends on the rate at which the driving field changes. 
We have demonstrated that in both universality classes 
the scaling exponents decrease logarithmically with in- 
creased driving rate AE/J. This is in sharp contrast 
to the cases where the dominant type of disorder is 
due to random vacancies in spin-systems Q and other 
avalanche-type dynamic systems driven at a finite rate 
p0{ , where the scaling exponents decrease and fractal di- 
mensions increase linearly with the driving rate. There- 
fore, occurrence of the logarithmic dependences of the 
scaling exponents on the driving frequency may serve as 
an additional test of the presence of the real random fields 
in ferroelectrics. 

Throughout this work we used the simulation algo- 
rithm at finite driving rates AE > 0, that mimics the 
realistic situation in experiments. The obtained set of 
scaling exponents, together with their dependences on 
the driving rate, belongs to one (of the two) universality 
classes, which implies a definite domain structure of the 
system, and hence the expected character of the switch- 
ing processes. Our results suggest that a crossover from 
one to the other universal behavior both in the switch- 
ing current noise and in the polarization reversal and 
relaxation can be achieved by varying the strength of 
pinning due to random electric fields. As mentioned 
above, in practice there are various ways to accomplished 
the crossover by reducing the pinning in ferroelectrics 
P, pJp7| . Here we have not discussed the properties of 
the dynamic phase transition separating these two scal- 
ing behaviors at a critical value / c , nor the critical expo- 
nents describing the depinning transition along the criti- 
cal field Hd(f) for / < f c . This questions require an ad- 
ditional study and proper finite size scaling analysis, for 
instance along the lines as done for the two-dimensional 
case in Ref. |28|l . Finally, by allowing a frequent ap- 
pearance of the unscreened depolarization fields a finite 
probability of negative pulses in the switching current 
builds up (b 7^ 0). The presence of negative pulses slows 
the switching process, interferes with the effects of finite 
driving rates, and alters the scaling properties of noise. 
A more complete analysis of the switching processes for 
non- vanishing probability b ^ is left for a future study. 

Finally, it should be stressed that, in spite of its sim- 
plicity, the proposed model seems reasonable for dynamic 
properties of three-dimensional uniaxial relaxor ferro- 
electrics. Whereas, in the single-crystalline normal ferro- 
electrics a more realistic modeling of the effects of inter- 
actions and depolarization fields on domain structure is 
required. Nevertheless, we believe that, when restricted 
only to the universal scaling behavior of the noisy com- 
ponent of the switching current, the presented results for 
the case of large domains may not depend significantly 



G 



on the details of interactions. Further research on the 
subject is necessary. 
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/ (pinning) 


Exponent 


3D 


2D 


/«/e 


T s 


1.78 


1.54 


(weak) 


Tt 


2.20 


1.83 


/ > fc 


T s 


1.66 


1.30 


(strong) 


Tt 


2.03 


1.47 



TABLE I. Scaling exponents for weak pinning (/ ~ f c ) 
and for strong pinning (/ > f c ) for 3-dimensional (3D) sys- 
tem at driving rate 4 x 10 -4 , results of this work, and for 
2-dimensional (2D) system at zero driving rate, from Ref. 
[25]. Estimated error bars are within ±0.03. 



(a) 
























/ / / 




3 -2 


-1 1 


2 ; 




Field 



0.8 0.9 1 1.1 
Field 




(d) 



I 



1200 



200 



800 



FIG. 1 

ing rate r • 



400 600 
Time 

Hysteresis loop for driv- 
0.01 and varying strength of disorder 



/ =4, 3.2, and 2 (inside out), (b) Sections of the first two 
loops from (a) for two different driving rates r =0.01 (dark) 
and r =0.08 (clear symbols), (c) Variation of the driving field 
with time and (d) the electrical Barkhausen noise (switching 
current) for finite probability of back flips b — 0.3 . Fig. 1 
from Ref. [18]. 



7 




FIG. 3. Cumulative distributions of switched domain size 
P(s, f) vs. size s (lower panel) and duration of switching 
P(t, f) vs. duration t (top panel) for varying strength of 
pinning / =2.2, 2.4, 2.8 and 3.2 (bottom to top, both pan- 
els). Driving rate is fixed AE/J = 0.001. Fitted values of 
three slopes, corresponding to full lines, are (bottom to top): 
t s - 1=1.029(4), 0.747(2), 0.568(5), on the lower panel, and 
T t - 1= 1.482(6), 1.145(7), and 0.94(2), on top panel. In total 
4.32 x 10 6 spins divided into 20 different samples was used 
in simulations. Distributions are normalized by the actual 
number of avalanches during the full cycle and the number of 
samples. Data are log-binned with ratio 1.1 . 
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FIG. 2. Part of the switching current signal j(t) against 
cumulative time t (lower panel) and the Fourier spectrum of 
the signal (top panel) for moderate pinning / = 3.6 . Slope 
of the fitted area <j> = 0.738 ± 0.007 . Only nonzero values of 
current are shown (see text). 
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FIG. 4. Cumulative distributions of switched domain size 
P(s, AE I J) vs. s (lower panel) and duration of switching 
P(t, AE I J) vs. t (top panel) for fixed pinning / = 3.2 and 
varying driving rate AE/J = 4 x 10~ 4 , 1 x 10~ 3 , 2 x 10~ 3 , 
4 x 10~ 3 and 8 x 10 -3 (bottom to top on both panels). Nor- 
malization to the actual number of avalanches results in the 
vertical shift of the curves. Scaling region indicated by the 
section of the straight line. Log-binning ratio 1.1 . 
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s 

FIG. 5. Same as Fig. 4 but for weak pinning regime 
/ = 2.4. Driving rates are (bottom to top) AE/J = 4 x 10~ 4 , 
1 x 1(T 3 , 2 x 1(T 3 , 4 x 1(T 3 , 8 x 1(T 3 and 1.2 x 1(T 2 . Fitted 
area on each curve is shown by solid line. 




Driving Rate AE/J 

FIG. 6. Scaling exponents for the distribution of switched 
area (squares) and duration of pulses (circles) against driv- 
ing rate AE/J. Filled and open symbols correspond 
to strong and weak pinning, respectively. Fit lines: 
Tx = A x - B x \n(AE/J), see text. 
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FIG. 7. Polarization reversal at field E/J — 3 (left part) 
and relaxation after field was set to zero (right part) against 
time for three values of pinning strength /, as indicated. 
Latter two lines were shifted to the right for better display. 
Fits according to P — 1 — C(t — a) - " (dotted lines), and 
P = A + Bexp[-(t/2 - 11) CT ] (full lines), see text. 
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